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We study the 2D Kondo insulators in a uniform magnetic field using quantum Monte Carlo simu- 
lations of the particle-hole symmetric Kondo lattice model and a mean field analysis of the Periodic 
Anderson model. We find that the field induces a transition to an insulating, antiferromagnetically 
ordered phase with staggered moment in the plane perpendicular to the field. For fields in excess of 
the quasi-particle gap, corresponding to a metal in a simple band picture of the periodic Anderson 
model, we find that the metallic phase is unstable towards the spin density wave type ordering for 
any finite value of the interaction strength. This can be understood as a consequence of the perfect 
nesting of the particle and hole Fermi surfaces that emerge as the field closes the gap. We propose 
a phase diagram and investigate the quasi-particle and charge excitations in the magnetic field. We 
find good agreement between the mean-field and quantum Monte Carlo results. 
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I. INTRODUCTION 

Kondo insulators, or heavy fermion semiconductors, 
are materials containing at least one atom per formula 
unit with a partially filled / or d shell and exhibiting 
properties similar to very narrow gap semiconductors. 
CeRhAs, CeRhSb, YB 12 , Ce 3 Bi 4 Pt 3 and SmB 6 are the 
most thoroughly investigated examples^. In the canon- 
ical model, the formation of the gap in Kondo insula- 
tors is a consequence of the hybridization between the 
conduction band and the effective /-electron level which 
gives rise to quasi-particle and spin-gaps at low temper- 
atures. Adopting a band picture one can close the gap 
by applying a high magnetic field, since the gap is on 
the meV scale. Although experiments on YB12 2 , SmB^ 
and Ce^Bi^Pt^ seem to support this simple picture, the 
exact nature of the field induced insulator to metal tran- 
sition as well as the role played by the strong correlations 
remains far from understood. 

Magnetic instabilities of the periodic Anderson model 
at half-filling have been studied extensively by slave bo- 
son mean field approximations^*^. A phase diagram as 
a function of the interaction strength was established and 
some thermodynamic and transport properties have been 
calculated. In these studies, the only effect of the mag- 
netic field is assumed to be the stabilization of the fer- 
romagnetically ordered state with respect to other mag- 
netic configurations. 

Carruzzo and Yu& studied the one dimensional, half- 
filled Kondo lattice in magnetic field using DMRG and 
bosonisation techniques. They found that although the 
spin gap closed at the critical field, the charge gap re- 
mained due to umklapp scattering They conclude that 
the ID half-filled Kondo lattice is insulating at all fields. 

The effect of the field in disordered Kondo insulators 
was treated by CPA in Ref HIol The authors find that, 
in the absence of magnetic ordering, the magnetic field 
induces the insulator to metal transition in the univer- 
sality class of density driven metal-insulator transitions. 
Based on scaling arguments the field dependence of the 



quasi-particle gap as well as the critical field as a function 
of temperature and impurity concentration were derived. 

In this paper we present a detailed study of the field 
induced quantum phase transition in 2D particle-hole 
symmetric models of Kondo insulators. We present a 
mean field calculation appropriate for the small-C/ limit 
of the periodic Anderson model. We find that the mag- 
netic field induces a phase transition from the paramag- 
netic insulator into a canted antiferromagnetic insulator 
which remains stable at all field strengths (until all the 
electrons in the system align with the field) . While zero- 
energy spin modes exist, we find that the field does not 
close the quasi-particle gap, if the lattice is bipartite, so 
that the metallic ground state is never induced by the 
field. We investigate more carefully the large- field limit 
using two effective models and reach essentially the same 
conclusion - on the bipartite lattice, the interaction is a 
relevant perturbation and the ground state remains in- 
sulating at all fields. The approximate treatments are 
complemented by a quantum Monte Carlo study of the 
particle hole symmetric Kondo lattice model in 2D. We 
find good agreement between the results. 

Recently, Beach and collaborators studied the effect 
of the magnetic field on the Kondo insulators using a 
large- N type mean field analysis of the Kondo lattice 
model and quantum Monte Carlo simulations 11 . They 
find that a large enough magnetic field induces a phase 
transition to a metallic ground state from the insulating 
canted antiferromagnetic state. The phase transition into 
the metallic state occurs when the / moments decouple 
from the conduction band, i.e. the hybridization mean 
field vanishes at a certain critical field. The question 
naturally arises whether this phase transition is real or 
possibly an artifact of the large- N mean field approach. 
Here we will show results which lead to a different conclu- 
sion: the insulating state induced by the magnetic field 
remains stable up to the full polarization of all electrons 
in the system, if the system is completely particle-hole 
symmetric. Thus the particle-hole symmetric Kondo in- 
sulator is an insulator at all fields. 
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The paper is organized as follows: In the next section 
we introduce the models used to describe Kondo insula- 
tors. In section IITT1 the phase diagram in the presence of 
the magnetic field is obtained using a mean field approx- 
imation for the half-filled periodic Anderson model. In 
section Hvl we present a discussion of the Kondo lattice 
model in high magnetic fields. In section the results 
of the Quantum Monte Carlo simulations are presented 
and compared with the mean field calculations. We sum- 
marize our results in section IVT1 and briefly comment on 
their relevance for the experimental systems. 

II. MODELS 

The canonical model used to describe the physics of 
the Kondo insulators is the periodic Anderson model 
(PAM)- 1 . The PAM Hamiltonian, including the uniform 
magnetic field in the z-direction is 

Hp AM = - ]T tijC^Cja+Cf flfiv+U n A n {i 
<i,j>,a icr i 

+ E( y ^ c ^ + H ' c -) - ^bB ■ + S?) (1) 

k,(T i 

Here all the symbols have their usual meaning. The PAM 
describes a two band system in which one band (conduc- 
tion electron, c band) is dispersive and uncorrelated and 
the other (/ band) dispersionless and strongly correlated. 
tij is the hopping matrix element in the c band and U 
the local Coulomb interaction in the / band. The two 
bands are mixed and the hybridization matrix element V 
controls the mixing strength. In the particle-hole sym- 
metric model that we consider in the following, tij = t 
for nearest neighbor sites on the square lattice and zero 
otherwise and e/ = —U/2. The magnetic field is coupled 
to the c and / electron spins only. The g factors of the c 
and / electrons are chosen to be the same, g c = gf = 2, 
for simplicity but choosing them differently would not 
change the qualitative aspects of our conclusions. In the 
following, the magnetic field is measured in the units of 
Zeeman energy. 

In the non-interacting (U — 0) case the ground state 
of the PAM is a paramagnetic band insulator with the 
quasi-particle gap A° p = ^(W/2) 2 + V 2 - W/2 ~ 
where W is the conduction electron bandwidth. In the 
field, the Zeeman splitting reduces the quasi particle gap. 
For fields larger than B c \ = A° p , the gap vanishes and 
th e ground stat e is metallic. In the fields beyond B c i = 
y/iW/2) 2 + V 2 + W/2 ~W + 2V 2 /W, all the spins are 
aligned with the field. The fully polarized ground state 
consists of two completely filled bands and is a trivial 
band insulator. 

Because of the particle-hole symmetry, the Fermi sur- 
faces of the spin up electrons and the spin down holes 



in the metallic state at intermediate fields are perfectly 
nested with respect to Q = (w, it). The staggered suscep- 
tibility in the plane perpendicular to the field diverges 
logarithmically as u> — > 0. This divergence makes the 
state unstable under perturbations coupling to the stag- 
gered magnetization. In particular, one expects that a 
staggered magnetization will be induced by any non-zero 
correlation on the / sites. The ensuing ordered state is 
a canted antiferromagnet, characterized by both m z and 
m x different from zero. 

When U is large enough (U/V ^> 1) to suppress 
charge fluctuations on the / sites, the low-energy physics 
of PAM is well described by the Kondo lattice model 
(KLM) 1 ^ 3 ., 

H KLM = -* E 4^a+jY,^-S{-2B z ^2(S^+Sr) 

(i.j)a i i 

(2) 

In the KLM, the charge fluctuations on the / sites are 
completely suppressed, / electrons are treated as spins 
and the hybridization is replaced by an antiferromag- 
netic exchange interaction between conduction electrons 
and / spins. Formally PAM and KLM can be re- 
lated by the Schrieffer- Wolff transformation 13 ' 14 , yielding 
J = 8V 2 /U. 

The zero-temperature, zero-field phase diagram of the 
2D particle-hole symmetric KLM has been well estab- 
lished by various numerical methodsi^iSiii. In the ab- 
sence of the magnetic field the ground state of the KLM is 
a paramagnetic insulator at large J/t. There is a quan- 
tum critical point at J/t ~ 1.4 and for small J/t the 
ground state is antiferromagnetically ordered. 

The large J/t paramagnetic state of the KLM is adi- 
abatically connected to the {7 = state of the PAM. In 
the particle-hole symmetric case, on finite lattices, this is 
guaranteed by theorems for the ground states of the two 
models^. 



III. MEAN FIELD ANALYSIS 

In this section we investigate the effect of the mag- 
netic field on the small-{7 PAM. In particular we want 
to investigate the spin density wave instability of the 
metallic state induced by the field in the non-interacting 
model. To this end, we perform the mean field decou- 
pling of the interaction term in Eq.Q by assuming the 
magnetization of the / spins to have a uniform compo- 
nent along the field axis and a staggered component in 
the plane perpendicular to the field, (S() — fhi with 
rfii = ((— Ym x , 0, m z ) . This yields the mean field Hamil- 
tonian (see the appendix lAl for details of the derivation), 
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Hmf = 5> fe - VoB)c\ a c k(y + J2(-P°)( B + Um z )flf ka + V^c^Aa + / fe W) 

- Um x ^2(ft +Q1 fki + fl+Qifc) + NU ( m l + ( 3 ) 



■^-[cos(fe x ) +cos(fc. y )]. The mean field Hamiltonian is quadratic in fermion operators 



with ^ = 1Q.) and e k — — g" i™°\"-xj -r ^\n y 

and is easily diagonalized by a unitary transformation. In the presence of the staggered magnetization, the Brillouin 
zone is halved and one finds 8 quasi particle bands; the particle bands 



££±(fc) = -j=[(B + Um z ) 2 + {Um x f + 2V 2 + {e k - Pa B) 2 ± 



± V((Um x ) 2 + (B + Um z f - (e k - p CT B) 2 ) 2 + 4V 2 [(e fc - 2 Pa B - p a Um z y + {Um x f 



1/2 



(4) 



and the hole bands related by, E% s (k) = —E° s (k). Note that the k dependence of the quasi-particle bands originates 
only from the dispersion of the conduction electrons. On a bipartite lattice, with e k+ Q = —t k , the quasi particle 
bands satisfy, E% a (e k ) = -E° s (e k+Q ). 

I 



In the ground state, the particle bands are empty 
and the hole bands are completely filled. To obtain the 
ground state energy, the expression 



E 9S =Y.Y;^)+ {Umx)2 t {Umz)2 (5) 



s.cr k 



u 



must be minimized with respect to m x and m Zl yielding 
the usual mean- field equations, 



8E, 



d(Um x 



0. 



(6) 



The prime on the summation sign in equation[3]indicates 
that the summation is to be taken over the magnetic 
Brillouin zone. 



A. Mean field phase diagram 

The minimization of the ground state energy for a 
range of U and B values was performed numerically and 
the obtained magnetization values are shown in Fig.^ In 
zero field, the system is paramagnetic at small values of 
U and antiferromagnetically ordered beyond U c ~ 1.25V. 
The staggered magnetization grows as (U — Uc) 1 ' 2 close 
to U c and tends to the fully saturated value m x = 1/2 as 
U -> 00. 

The magnetic field applied to the system reduces the 
value of U at which the magnetic instability occurs. The 
phase boundary can be obtained by solving 



d 2 E 



U C (B) 



d{Um x f 



(7) 



where E = ^ s=± a=u J2k E h,si^)- At small fields, 
B <C B c i, the critical interaction strength falls off as the 
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FIG. 1: Staggered (left) and parallel magnetizations vs. B 
and U for W = 1, V = 1 obtained by numerically minimizing 
the mean field equations. The grayed out plateau marks the 
paramagnetic phase. 



square root of the field, U C (B) — U c (0) oc —yB, as ex- 
pected in the mean field approach. At the phase bound- 
ary one finds the usual mean field critical exponents for 
the staggered magnetization, oc (U — U^B)) 1 / 2 and 
m x s oc (B — B^U)) 1 / 2 . After the initial rise, m x goes 
through a maximum and falls of exponentially in large 
fields. On the phase boundary, the parallel susceptibility 
vanishes (it is zero in the paramagnetic phase, since the 
spin excitations are gaped) . Close to the phase boundary, 
it behaves like m z oc (B — B c (U)) a , with a > 1. 

The right hand side of the Eq.Q is proportional to 
the static staggered susceptibility of the Um x = state, 
Xo (Q)- This can be expressed using the familiar Lind- 
hardt formula which in the case considered here reduces 
to 



x«r(Q) = -E 



,f{El_{k))-f{El_{k)) 



El_(k)-El_(k) 



(8) 



U m x —0 



In small fields, the quasi-particle gap provides a cut-off 
for the denominator in the sum on the right hand side 
and xt~(Q) IS nn he. When the field closes the gap, 
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FIG. 2: Quasi-particle gap (dots), staggered magnetization 
(dashed line) and total (/ + c) parallel magnetization (thin 
solid line) vs. magnetic field for W = 4, V = 1, U = 1. In the 
paramagnetic phase the gap decreases linearly with field. At 
large fields the gap follows the staggered magnetization. 



the denominator vanishes along the Fermi surface (Fermi 
lines), determined by the equation, 



e k = ±e = ±- 



V 2 - BUm 2 



B 



B + Um z 



(9) 



Consequently, the staggered susceptibility diverges loga- 
rithmically in the field B > A^ p and there is no finite U c . 
The system is ordered for any finite interaction strength. 
The divergence of \o (Q) is a direct consequence of the 
perfect nesting, _(k) = —E^ _(k), and is found at all 
fields, if the conduction electron hopping is constrained 
to a bipartite lattice and the system is half filled. 18 

The behavior of the staggered magnetization at small 
U can be found by solving the mean field equations to 
leading logarithmic order in U m x . The details of the cal- 
culation are described in the appendix lAl The resulting 
expression for the magnetization is 



m x oc exp 



(B + Um, z f + V 2 {B + Um z ) 2 



V 4 P0 U 



(10) 



where po is taken to be the density of states on the 
m x = Fermi surface. It is interesting to note that 
the expression ((T0"|> is valid for the large field region and 
for the large U region with B 3> 8V 2 /U. 

B. Quasi-particle spectrum 

The field dependence of the quasi-particle gap for a 
fixed value of the interaction is shown in the Fig. [3 In 
the paramagnetic phase, the gap decreases linearly with 
the field. In the ordered phase, the quasi-particle gap 
is proportional to the staggered moment and follows the 
same exponential dependence for large fields. It is impor- 
tant to realize that the quasi-particle gap always remains 
finite, so that the system is insulating. 
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FIG. 3: Poles of the a =\ conduction electron spectral func- 
tions along the high symmetry lines of the Brillouin zone for 
field and interaction strengths indicated in the plots. The 
width of the lines indicates the weight in the pole. The plots 
on the right show the contours of constant quasi-particle gap 
magnitude with dashed lines indicating the position of the 
Fermi surfaces of the m x — state. 



The spectral functions for the electrons in the mean 
field model show infinitely sharp peaks at the quasi- 
particle band energies. The poles of the a —1 electron 
spectral function in the ordered phase at various values 
of the field and interaction strengths are shown in Fig. [3] 
The width of the lines in the figure indicate the weights in 
the corresponding poles. Note that the gap at the Fermi 
surface is always finite, even though the exponentially 
small scale is not immediately apparent in the plots. The 
contour plots show the quasi-particle gap size in the Bril- 
louin zone. The gap minima are indicated by the dashed 
lines in the contour plots. The location of the gap min- 
ima indicates also the position of the Fermi surfaces of 
the m x = 0, metallic state. 

It is interesting to observe the change in the character 
of the quasi-particles at the gap minima as the field and 
the interaction strength are varied. At small U and B ~ 
A° p , the low-energy quasi-particles are "heavy" and the 
minimum of the gap lies near the zone center. As the 
field is increased, the minimum moves towards the zone 
diagonal and the quasi-particles become more and more 
c like. 

When the gap minimum reaches the zone diagonal, the 
magnetization of the system along the field direction is 
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exactly one half of the fully saturated value. The field 
strength at which this happens, B 1 / 2 , depends on the 
interaction strength and can be obtained by setting eo = 
with the limiting behavior: 



(11) 



In the large U limit, Bi/ 2 sets the energy scale at which 
the / electrons align with the field. 



IV. LARGE B LIMIT OF THE KONDO MODEL 

For large values of U, in the fields B > Bi/ 2 , the / 
electrons are almost completely aligned with the field. 
At the mean field level, the poles of the spectral function 
corresponding to the charge fluctuations on the / sites 
move towards ±(B + Um z ), i.e. far from the Fermi level. 
In the large- N mean field theory this eventually results 
in the complete decoupling of the / electrons from the c 
band and the decoupled metallic state obtainsiii 

In the particle hole symmetric case, the Fermi surfaces 
of the metallic state are perfectly nested. The perfect 
nesting makes the metallic state unstable at all fields in 
the small U limit of the PAM. We will now demonstrate 
that also in the limit of large Coulomb interaction, i.e. 
for the KLM, the same instability arises. 



A. Effective Hamiltonian approach 

We consider the KLM Hamiltonian in the large mag- 
netic field B J. In the magnetic field, the J = 
ground state of the KLM is non degenerate and is given 

by IV;) = IW T 4 T Uk<k Fl Ai IL 4 1)- Flipping the / 
spin is an excitation with a gap given by the Zeeman en- 
ergy. A canonical transformation approach can be used 
to generate an expansion in (J/ B) around the J = 
ground state. The effective Hamiltonian governing the 
low energy dynamics of the system is given by (the de- 
tails of the derivation in the appendix IE|) . 



H 



E 



(ejt - p a B)cl c k 



(12) 



where U — j^, e k is the dispersion of the original conduc- 
tion electron band and B = (B — J/4 — U/2) is the effec- 
tive magnetic field. In this effective model, the spin flip 
interaction between the conduction band and the fully 
polarized /-spin background of the KLM, has been re- 
placed by a contact interaction between the c electrons 
and the / spins have decoupled from the dynamics. 

If the conduction electron band is particle hole sym- 
metric, so that ek+Q = — £fc, the spin up hole and 
the spin down electron Fermi surfaces of the effective 
model are perfectly nested. Any non zero U there- 
fore induces magnetic ordering in the plane perpendic- 
ular to the applied field. A mean-field decoupling, with 



(s*j) = ((— ) % m x , 0, m z ), analogous to the one performed 
in section ITTTI yields the quasi particle bands 



E a , ± (k) = ±J(e k - PlJ (B + Um z )) 2 + (Um x ) 2 (13) 



and the mean-field equation determining m x 



2 

U 



p(e)de 



1 



1 



£ T ,+ (e) E L+ (e) 



(14) 



In high magnetic fields, the up and down spin Fermi sur- 
faces are well approximated by circles of radii W — B 
centered at (tt, tt) and (0,0), respectively. We therefore 



can set p(e) 



Po 



~ ? to obtain 



Um x oc 2(W — B — Um z ) exp 



PoUj 



(15) 



The staggered magnetization and the quasi-particle gap 
are finite for any finite U, as long as B + Um z < W. It is 
easy to see that, B + U m z — W is just the condition for 
system to fully polarize. This means that the staggered 
magnetization vanishes only in the completely polarized 
system. The completely polarized phase is a trivial in- 
sulator. The metallic state is, therefore, never obtained 
in the particle-hole symmetric case and is a bad starting 
point for the perturbation expansion. 



B. Classical spins mean field 

We have seen that a large magnetic field suppresses 
both charge and spin fluctuations on the / sites. The 
physics of the high-field phase will, therefore, be well 
described by the KLM in which the / spins are replaced 
by an array of statically arranged classical spins. Let the 
spin configuration be 



S, 



sin 9 cos Qri 
— sin 9 sin Qr^ 
cos 9 



(16) 



With Q — (tt, tt), this corresponds to the same choice of 
/ magnetization as in section fllll with 1/2 sin 9 = m x 
and 1/2 cos 9 = m z , so that the system is fully polarized 
for 9 = 0. The problem is now reduced to one of the non- 
interacting conduction electrons in an external magnetic 
field, described by the Hamiltonian, 

H = ( £ fe - P*B) c L c ka - NB cos 9 

Jsin9 s-^ ( t 4. \ 

+ 1^ (4-QT c H + 4+Qi c fcT J . (17) 

k 

with B = B — J/4 cos 9. This is easily diagonalized to 
find the quasi-particle bands (k in the magnetic Brillouin 
zone) , 



E„, ± {k) = ±J{e k - Pa B)2 + (J/4) 2 sin 2 9. (18) 
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In the ground state the "-" bands are completely filled. 
Minimizing the ground state energy and assuming the 
same circular Fermi surface approximation as in the pre- 
vious subsection, the mean field equation determining the 
angle 9 is obtained as, 



B tan 9 = 



Pa 



-w ' 



E 



l\ s,„f 



-p a (e - p a B)^ ta.nl 



E, 



(19) 



For sin 9 <C 1 (/ moments almost aligned with the field) 
we obtain, 



8B Ap B + 8 



Jm x oc y (W - B)(W + B) exp 



PoJ 2 



PoJ 



(20) 

The staggered magnetization vanishes only when B = W. 
It is easy to see that this is exactly the condition for the 
system to fully polarize. The full polarization field is 
equal to W — J/4 and agrees with the one obtained us- 
ing the effective Hamiltonian. The dominant, small J 
exponential dependence of the staggered magnetization 
m x oc exp[— 8B/(p J 2 )] also agrees with the one obtained 
in the previous section. The subleading 1/(pqJ) correc- 
tion to the exponent arises because the Zeeman energy 
of the / electrons has now been taken into account. As 
the quasi-particle gap is proportional to the staggered 
magnetization, the system stays insulating at all fields. 



V. QUANTUM MONTE CARLO 

In this section we present QMC simulations of the 
Kondo lattice model in the magnetic field. As in the 
zero field case, the sign problem may be avoided only 
for particle-hole symmetric conduction bands. To com- 
pare with the mean-field results we adopt a projective 
approach in which the ground state, |^o)i 1S filtered out 
from a trial wave function, \^t), satisfying (^oI^t) 0. 
In this algorithm, the ground state expectation value of 
an observable O is estimated via: 



(*o|0|* > ,. (y T \e- eH/2 Oe~ 0H / 2 \y T ) 

— : ; = lim — 

(*o*o) 



(21) 



In the QMC, we evaluate the right hand side of the above 
expression at finite values of O and then extrapolate to 
infinite values. The details of the algorithm - in particu- 
lar the sign free formulation - has been described exten- 
sively in Ref. ITrj Since we are working in the canonical 
ensemble, the total magnetization 



ATf + Nj - A j - N l f 

K 



(22) 



with N u the number of unit cells, is fixed during the 
simulations. 



By measuring time displaced correlation functions, we 
can extract quasi-particles as well as spin gaps. Consider 

(*o|5-(-g,r)5+(?,0)|* > = 
5^Kn|S+(g)|*o>| 2 e- T ( B »(«' JV ' s * +1 >- Bo ^^) (23) 

n 

where E n (q,N,S z ) are eigenstates of H with momen- 
tum q, particle number N and total z-component of 
spin S z . From the large r behavior of the above cor- 
relation functions, we can extract the energy difference 
Eo(q, N, S z +1)—Eq(N, S z ) from which we can determine 
the spin gap: 

A sp (q) = E (q, N, S z + 1) - E (N, S z ) - h (24) 

where h = [E {N, S z + l)- E (N, S z - 1)] /2. In the 
same manner, we compute the quasi-particle gap from 
the single-particle imaginary time displaced Green func- 
tion. 

A qp (k) = E (k, N + 1, S z ) - E (N, S z ) - y. (25) 

with chemical potential: \x = [E (N + 1, S z ) — E (N — 
1, S z )]/2. In Fig.0|we plot the gaps as a function of total 
magnetization at J ft = 2.0. 

In the zero field case, the Kondo insulating state with 
finite quasi-particle and spin gaps is realized. At fi- 
nite magnetizations and according to the mean-field ap- 
proach, we expect a canted antiferromagnetic state and 
hence no spin gap. Fig.0]confirms this point of view: the 
spin gap drops to zero at all finite values of the magne- 
tization within the accuracy of the numerical simulation 
and the equal time spin-spin correlations in the plane 
perpendicular to the magnetic field show long range an- 
tiferromagnetic order. In the mean-field approach, mag- 
netism stems from a Fermi surface instability and due to 
perfect nesting opens a gap on all the Fermi line. The 
QMC results of Fig. 0] are consistent with this predic- 
tion, since, as apparent, the quasi-particle gap survives 
at finite magnetizations. 

To further compare the mean-field approach to the 
QMC we have used the Maximum Entropy method to 
obtain the single particle spectral function A(k,u>). At 
zero field in the Kondo insulating phases, the dominant 
features of the spectral function are well described by 
hybridized bands (solid vertical lines in Fig. |SJi). From 
the single particle occupation number a — (ct C£ CT ) 

listed on the left hand side of Figs [SJl-c we can extract 
the total weight under the "photoemission" (uj < 0) and 
the "inverse photoemission" (lo < 0) spectra since 

/•OO 

dujA a (k,oj) =im Sa / dujA a (k,uj) =7r(l-n jEtr ). 

_ j 'Jo 

In particular, one sees that the photoemission (inverse 
photo-emission) spectrum in the vicinity of k = (tt, 7t) 
(k = (0, 0)) has a small weight. Those heavy bands stem 
from the Kondo screening. 
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FIG. 4: (a) Quasi-particle and spin gaps as a function of the magnetization, (b) x-component of the f-spin-spin correlations at 
the largest distance, (L/2, L/2), on the L x L lattice. The Fourier transform of this quantity, 2~2 f! e''^ r {m£(0)ml(r)) , is peaked 
at the antiferromagnetic wave vector, Q = (tt, 7r). (c) Size scaling of the quasi-particle gap at M z — 0.25,0.5. The data is 
consistent with a finite value of A qp in the thermodynamic limit, (d) Size scaling of S£ (Q = (n, n)) = '}2 i ,e l( ^' r (ml(f)ml(0)} for 
M z = 0.5, 0.25 As apparent, the data is consistent with S[(Q) ~ L 2 thus signaling the presence of long-range antiferromagnetic 
order perpendicular to the applied field. 



In the mean field approach, the effect of a magnetic 
field is to shift the spin down band up in energy until 
it ultimately crosses the Fermi surface, thus generating 
a metallic state. This metallic state is however unstable 
due to the underlying particle-hole symmetry. In Fig. 
we compare the finite field results with a rigid shift of the 
hybridized bands. The data are compatible with the in- 
terpretation that the down spin band has indeed crossed 
the Fermi surface but that at the crossing point the mag- 
netic instability opens a gap. Furthermore on the photo- 
emission side around the (tt, tt) point we see a very weak 
feature which we can identify as a shadow of the up band 
which has dropped below the Fermi energy in the vicinity 
of the k = (0, 0) point. 

Breaking of the spin symmetry by the magnetic field 
suppresses Kondo screening. We hence expect the weight 
of the features in the spectral function stemming from 
Kondo screening to be suppressed as a function of grow- 
ing magnetization. For example, consider the inverse 
photoemission at k = (0, 0) in Fig. |SJ As apparent the 
weight of this feature is reduced as a function of growing 
values of M z and will vanish when the /-spins become 
fully polarized. In our simulations, the /-spins are never 
completely aligned with the field as long as the system is 
not completely polarized, i.e. as long as Nj + N. < N u . 



This supports the conclusion that the metallic state is 
never induced by the field. 



VI. CONCLUSION 

We studied the magnetic field induced quantum phase 
transition in the 2D particle-hole symmetric Kondo insu- 
lators using: i) a mean field approximation appropriate 
in the small U limit of the periodic Anderson model, ii) 
two mean field approximations appropriate in the large 
field limit of the Kondo lattice model and iii) a quantum 
Monte Carlo simulation of the particle-hole symmetric 
2D Kondo lattice model in the field. 

We find a magnetic field induced quantum phase tran- 
sition from a paramagnetic insulator into a canted anti- 
ferromagnetic insulator ground state. In the particle-hole 
symmetric case we studied, the antiferromagnetism can 
be understood as a spin density wave type instability of 
the perfectly nested quasi-particle Fermi surfaces that 
would arise in the field in the absence of interactions. 
Because of the perfect nesting, any finite interaction is a 
relevant perturbation and results in a finite quasi-particle 
gap. Consequently, the ground state of the interacting 
system remains insulating in all fields. 
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FIG. 5: Single particle spectral function in the down-spin sector as a function of magnetization as obtained from analytical 
continuation of the QMC results with the Maximum Entropy method. The solid vertical lines correspond to the pole position 
of hybridized bands (see text). On the left hand side of each Fig. we have listed the single particle occupation number . 
from which one can determine the weight under the spectral function (see Eq. 1261 . 



We conclude that the recently proposed insulator to 
metal transition induced by the fieldii is likely to be 
an artifact of the large- N approximation to the Kondo 
lattice model in the particle-hole symmetric case. If, 
however, the particlc-holc symmetry is violated a field- 
induced metal-insulator transition is possible in certain 
parameter ranges. 

We find that the qualitative features of the phase dia- 
gram as well as of the quasi-particle excitations are well 
described by a simple mean field approximation to the 
periodic Anderson model. The magnetic field explicitly 
breaks the spin rotation symmetry and suppresses the 
charge fluctuations on / electrons, essentially by fully 
polarizing the / band. 

The band structure of the Kondo insulators is non- 
trivial and deviations from particle-hole symmetry are 
to be expected in the real materials. In the absence of 
perfect nesting there would be a critical value of the field, 
controlled essentially by the nesting mismatch, at which 
the gap will close on some parts of the Fermi surface. 
Therefore one expects to eventually find a metallic state 
induced by the field. Finally we would like to mention 
that the conclusions we draw here are valid also for the 
three-dimensional systems, where the same kind of nest- 
ing features would appear for perfect particle-hole sym- 
metry. 
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APPENDIX A: MEAN FIELD DECOUPLING OF 
THE HUBBARD TERM IN THE CANTED 
ANTI-FERROMAGNETIC PHASE 



We want to decouple the interaction term in the PAM, 
in the presence of a canted staggered magnetization. To 
this end, we select the spin quantization axis at each 
site to point in the direction of the local magnetiza- 
tion, fhi. Using the operator identity f^JiufiiJiU = 



l/4(4/i Ti + fljiU? ~ WflufiU ~ fli 
Ti (li) denotes the spin with respect to the local quanti- 



ze Au) 2 , wh ere 
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zation axis, the interaction term can be decoupled as 



N- 



Un 2 



U J2 ~fliM) + u\M\ 2 

i 

+ Fluct. (Al) 

where rif = + fjifil) is the average occupancy of 

the / site and "Fluct." denotes the terms neglected in the 
mean-field approximation. After the decoupling, a spin 
axis rotation to a common quantization axis (given by the 
direction of the external field) is performed using R = 
cxp —i/^i ' S(\ with Oi being the vector pointing 



along B x m; and of magnitude equal to the angle between 
B and fhi. Since, 

1^1^(4/^ - f} u hu)R ] =m-S{, (A2) 

this yields the mean field Hamiltonian, Eq. [31 of sec- 
tion [fa] 

To obtain the behavior of the staggered magnetization 
in fields B > A® p and at small U we need to solve the 
mean field equations [S] in the limit when Um x — » 0. As 
the only k dependence of the quasi-particle bands comes 
through the k dependence of the conduction electron en- 
ergy, the summations over k are readily transformed into 
integrals over the conduction electron energy, thus yield- 
ing (p(e) is the conduction electron DOS), 



d 



d(Um x ) 



-w 



p(e)de(El 



<±(e)) 



Urn, 
> 

' u 



■ = 0. 

For U m x — > 0, the dominant contribution to the integral 
comes from the band crossing the Fermi surface at eo 
(given by Eq. 0. For A° qp < B < V, this is E^ + . For 
small U and close to to we can write, 



with 



r 2 



a 



(B + Um z ) 2 + V 2 

(B + Um z ) 2 
(B + Um z ) 2 + V 2 ' 



(A4) 

(A5) 
(A6) 



The logarithmically divergent part of the mean field 
equation can now be written as (neglecting the non- 
divergent contributions), 



1 

U 



a 



de'p(e + e') 



( aUrn x \ 

V f> J 



(A7) 



where we have introduced a cutoff e c which does not in- 
fluence the exponential dependence. Eq. I1UI in the text 
now follows by elementary integration. 



APPENDIX B: EFFECTIVE HAMILTONIAN IN 
THE LARGE FIELD 

The ground state of the J = KLM in the magnetic 
field is non-degenerate. All the / spins are polarized in 
the direction of the field. The / spin flip is an excitation 
with an energy gap of the size of the Zeeman energy. Let 
P n denote the portion of the Hilbcrt space with n lo- 
calized spins pointing opposite to the magnetic field and 
let V n be the corresponding projector. In large magnetic 
fields, B ^> J, one expects the ground state and the low 
lying excitations to lie dominantly in Pq and have only 
small components in the P n , with n > 0, subspaces. We 
split the Hklm into a P n diagonal part 

Ho = $> fcff -p a B)c\ a c ka - 2Bj2p*Si 



J 



J2S?(4^-4 lCil ) (Bl) 



and the spin flip part 



Let S be a Hermitian operator, such that 
[H ,S] = -V. 



(B2) 



(B3) 



The effective Hamiltonian, obtained by applying the 
canonical transformation 



H = e b H KLM e'- 



(B4) 



has no matrix elements between the states in Pq and 
Pn>o of order less then 0(J(J/B) 2 ). We can therefore 
obtain the low-energy dynamics of the original problem 
correctly to order J 2 / B) by considering only the Pq part 
of the Hilbert space and the Hamiltonian, 



H F " =V e b H KLM e- b P Q . 



(B5) 



By expanding and rearranging the exponentials one ob- 
tains 



H^°=Vo ( H + Uv,S] 



To 



(B6) 



Using Eq. (|B3(I it is easy to obtain the matrix elements of 
S between the eigenstates of Ho from which the operator 
form easily follows, 

S = AB ^ ( 5 « +C Ul CfcT ^ S 7 c Ut c h) > (B7) 



kq 



thus yielding, 

j 2 \ 



[V,S] = 4BN K c feT Cfc -^' 5 M"i 



(B8) 



qkk' 
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Evaluating the commutator and projecting on the Pq sub- 
space, yields 

j2 ^ 

V [V,S]Vq = -T^T7^ C !'i Cfe '-<?T c !T C fc-9l- ( B9 ) 



kk' q 



The last term in the above equation is a constant and can 
be dropped, the sf term is the contribution to the effec- 
tive uniform magnetic field. This completes the deriva- 
tion of the effective model described in section Hv] 



Substituting into Eq. i|B6(l and neglecting the terms cor- 
responding to the dots, which are of order 0(J 3 /B 2 ) one 
obtains the effective Hamiltonian, 



72 j2 j2 w 

H = Ho + ^J2st+ — ^n i ,n il -—. (BIO) 
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